Improving efficiency of supercontinuum generation in photonic crystal fibers by direct 

degenerate four- wave- mixing 
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We numerically study supercontinuum (SC) generation in photonic crystal fibers pumped with 
low-power 30-ps pulses close to the zero dispersion wavelength 647nm. We show how the efficiency is 
significantly improved by designing the dispersion to allow widely separated spectral lines generated 
by degenerate four-wave-mixing (FWM) directly from the pump to broaden and merge. By proper 
modification of the dispersion profile the generation of additional FWM Stokes and anti-Stokes lines 
results in efficient generation of an 800nm wide SC. Simulations show that the predicted efficient 
SC generation is more robust and can survive fiber imperfections modelled as random fluctuations 
of the dispersion coefficients along the fiber length. 



I. INTRODUCTION 

After the first observation of a 200-THz supercontin- 
uum (SC) spectrum of light in bulk glass [1,2], much 
has been done on the understanding and control of this 
process [3] . A target of numerous experimental and the- 
oretical investigations has been the improvement of the 
characteristics and simplification of the technical require- 
ments for the generation of a SC [3]. The first experi- 
ments in bulk materials, based on self-phase modulation 
(SPM), required extremely high peak powers (> 10MW). 

New techniques based on the use of optical fibers as a 
nonlinear medium for SC generation allowed lower peak 
powers to be used due to the long interaction length and 
high effective nonlinearity [4-7]. However, the necessity 
to operate near the wavelength for zero group velocity 
dispersion, restricted the SC generation to the spectral 
region around and above 1.3/xm. The use of dispersion- 
flattened or dispersion-decreasing fibers as nonlinear me- 
dia for SC generation resulted in a flat SC spanning 1400- 
1700nm [8,9] and 1450-1650nm [10], respectively. The 
spectrum was still far from the visible wavelengths and 
in some cases very sensitive to noise in the input [10]. 

Photonic crystal fibers (PCF) and tapered fibers over- 
come these limitations. The unusual dispersion prop- 
erties and enlarged effective nonlinearities make them 
a promising tool for effective SC generation [6]. PCFs 
and tapered fibers have similar dispersion and nonlin- 
earity characteristics and they have the advantage that 
their dispersion may be significantly modified by a proper 
design of the cladding structure [11-13], or by chang- 
ing the degree of tapering [14] , respectively. Using kilo- 
watt peakpower femtosecond pulses a SC spanning 400- 
1500nm has been generated in a PCF [14] and in a ta- 
pered fiber [15]. The broad SC was later explained to be 
a result of SPM and direct degenerate four-wave-mixing 
(FWM) [16]. 

However, high power femtosecond lasers are not nec- 



essary, - SC generation may be achieved with low-power 
picosecond [6,7] and even nanosecond [17] pulses. Thus 
Coen et al. generated a broad SC in a PCF using picosec- 
ond pulses with sub- kilowatt peakpower and showed that 
the primary mechanism was the combined effect of stim- 
ulated Raman scattering (SRS) and parametric FWM, 
allowing the Raman shifted components to interact effi- 
ciently with the pump [6]. 

Using 200 fs high power pulses and a 1cm long tapered 
fiber, Gusakov has shown that direct degenerate FWM 
can lead to ultrawide spectral broadening and pulse com- 
pression [16]. We consider how the direct degenerate 
FWM can significantly improve the efficiency of SC gen- 
eration with sub-kilowatt picosecond pulses in PCFs, and 
go one step further by optimizing the effect through engi- 
neering of the dispersion properties of the PCF. We show 
that by a proper design of the dispersion profile the di- 
rect degenerate FWM Stokes and anti-Stokes lines can 
be shifted closer to the pump, thereby allowing them to 
broaden and merge with the pump to form an ultrabroad 
SC. This significantly improves the efficiency of the SC 
generation, since the power in the Stokes and anti-Stokes 
lines no longer is lost. In particular we optimize the SC 
bandwidth by designing the dispersion profile to generate 
additional Stokes and anti-Stokes lines around which the 
SC spectrum can broaden. 

External perturbations and different types of imper- 
fections lead to fluctuations of the fiber parameters along 
the length of the fiber. Fluctuations in fiber birefringence 
[18,19], dispersion [21,22], and nonlinearity [23,24] has 
been investigated to understand their influence on differ- 
ent regimes of light propagation. As parametric processes 
require phase matching, the effectiveness of the FWM 
could be strongly influenced by random fluctuations of 
the dispersion. Indeed Coen et.al. in their PCF exper- 
iments with low-pump picosecond pulses at 647nm [6,7] 
and nanosecond pulses at 532nm [17] explains the ab- 
sence of frequencies generated by direct degenerate FWM 
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from the pump, by the large frequency shift and the vi- 
olation of the required phase-matching condition due to 
fiber irregularities. We analyze the influence of a ran- 
dom change of the dispersion coefficients along the fiber 
on the process of SC generation and find that the gener- 
ation and merging of the direct degenerate FWM Stokes 
and anti-Stokes waves with the pump could be robust 
enough to survive fiber imperfections, and thus a signifi- 
cant improvement of the process of SC generation should 
indeed be possible in real PCFs. 



II. THEORETICAL MODEL AND FIBER DATA. 

We study numerically the SC generation process us- 
ing the well known coupled nonlinear Schrodingcr (NLS) 
equations that describe the evolution of the x- and y- 
polarization components of the field for pulses with a 
spectral width of up to 1/3 of the pump frequency [6,25], 



The phase-mismatch for degenerate FWM of two pho- 
tons at the pump frequency is: A/3 = (3 S + (3 as — 2(3 P + 
27(1 — fn)Ip [27], where I p is the peak power. In the 
frequency domain we have: 

A/3 = 2(tt 2 § + fi 4 § + + 7 (1 - f R )I p ), (3) 

where f2 = u p — uo s = u as — u p . The gain g of the direct 
degenerate FWM [7,27] is: 
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Here the complex fields Aj = Aj(t,z) with j = 1,2 are 
given by A\ = E x and A2 = E y exp(iS(3z), where E x 
and E y are the envelopes of the real linearly polarized 
x- and y-components. The time r = t — z/v is in a 
reference frame moving with the average group velocity 
v^ 1 = (?; 2 T 1 + v~ 1 )/2, z is the propagation coordinate 
along the fiber, fi is the fiber loss, 5(3 — (3 X — (3 y = uj^brij c 
is the phase mismatch due to birefringence 5n = n x — n y , 
and A = (uj 1 —Vy 1 ) is the group velocity mismatch be- 
tween the two polarization axes. The propagation con- 
stant (3{lo) is expanded to 8 th order around the pump 
frequency uj p with coefficients (3k keeping P2-7 same for 
x- and y-linearly polarized components, 7 is the effec- 
tive nonlinearity, Jr is the fractional contribution of the 
Raman effect, and finally * denotes complex conjugation. 

This model accounts for SPM, cross-phase-modulation, 
FWM, and SRS. For the Raman susceptibility Hr we 
include only the parallel component, as the orthogonal 
component is generally negligible in most of the fre- 
quency regime that we consider [26]. The Raman sus- 
ceptibility is approximated by the expression [27]: 
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where 17 = 12.2fs and T2 = 32fs. Furthermore, f R = 0.18 
is estimated from the known numerical value of the peak 
Raman gain [27]. 
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where I is the power of the frequency component around 
which the degenerate FWM process takes place. 

We consider the same PCF and use the same nu- 
merical and experimental data as in [6], kindly pro- 
vided by S. Coen. We pump along the slow axis with 
30ps sech-shaped pulses of I p = 400W peak power and 
pump wavelength X p = 647nm. The PCF has core area 
A CO re = 1.94,um 2 , birefringence 5n = 1.9 • 1CT 6 , and non- 
linearity 77,2 = 3 • 10~ 20 m 2 /W. We consider six different 
dispersion profiles, for which P2-7 are given in Table I, 
where case dl corresponds to the PCF used in [6]. Note 
that cases d4+d6 and d5 have two and three sets of Stokes 
and anti-Stokes waves, respectively. 

Here we are interested in studying the separate effect of 
different dispersion profiles (i.e. different values of fii-i)- 
We therefore keep all other fiber and pulse parameters 
constant and neglect the frequency dependence of the 
effective area A e ff and the loss /1. Thus we assume a 
uniform loss of /j, = O.ldB/m and approximate A e ff with 
the core area, so that 7 = 2n 2 / (X p A core ) = 0.15(Wm) _1 . 

TABLE I. Dispersion coefficients (32 [ps 2 /km], 
(3a [10 _5 ps 4 /km] and (3q [10 _10 ps 6 /km] for dispersion profiles 
dl-d6, with corresponding dispersion at the pump wavelength 
D(X P ) [ps/nm • km], zero dispersion wavelength X z [nm] and 
Stokes A s [nm] and anti-Stokes A as [nm] wavelengths. Fixed 



coefficients: (3z = 
and fi, = 1.2 • 10" 



5.1 • 10~ 2 ps 3 /km, p B = 1.2 ■ 10- 7 ps 5 /km 
3 ps 7 /km 
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fa 
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(3 6 


A z 


D(X P ) 


A s 




dl 


7.0 


-4.9 


-1.8 


677 


-31.6 


1110 


457 


d2 


14 


-34.4 


-0.04 


697 


-62.3 


852 


521 


d3 


1.0 


-2.5 


-3.3 


652 


-4.5 


852 


521 


d4 


-0.28 


0.05 


0.29 


647 


1.3 


1062 

852 


465 
521 


d5 


-1.01 


2.14 


-2.84 


643 


4.54 


1104 
894 
751 


458 
507 
569 


d6 


-1.3 


-2.6 


58.8 


641 


5.9 


800 
720 


543 
587 



2 



For our numerical simulation, we use the standard sec- 
ond order split-step Fourier method, solving the nonlin- 
ear part with a fourth order Runge-Kutta method us- 
ing a Fourier transform forth and back and the convolu- 
tion theorem. Except where otherwise stated, we use 2 17 
points in a time window of T = 236ps, giving the wave- 
length window (405 — 1613nm). The propagation step is 
Az = 43/im. In our longest simulation out to L = 3.7m 
the photon number is conserved to within 5% of its ini- 
tial value. An initial random phase noise seeding of one 
photon per mode is included as in [6]. All the presented 
spectra have been smoothed over 32 points. 



III. NUMERICAL ANALYSIS. 

We first simulate SC generation using the same fiber 
as in [6], i.e., using the dispersion profile dl. Due to 
our large spectral window (405 — 1613nm), we see in 
Fig. 1 (left) the emergence of direct-degenerate FWM 
Stokes and anti-Stokes waves at the predicted wave- 
lengths A s = llOOnm and A as = 458nm, for which the 
phase matching condition (3) is satisfied. From the stan- 
dard expressions given in [27] we find the maximum di- 
rect degenerate FWM gain (g) to be twice the maximum 
SRS gain, which explains why the FWM Stokes and anti- 
Stokes components appear before the SRS components. 
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FIG. 1. Dispersion profiles dl (left) and d2 (right). Top 
row: phase mismatch A/3 for direct degenerate FWM (dashed 
line) and spectrum at L=10cm (solid line). Bottom row: spec- 
trum at L=20cm, 30cm, lm, 2m, and 3.7m (down to up). 

For a given peak power, the loss and temporal walk- 
off of the PCF gives the maximum distance L max over 
which nonlinear processes, and thus the SC generation 
process, are efficient. From Fig. 1 (left) we see that 



after the FWM Stokes and anti-Stokes components are 
generated they broaden in the same way as the central 
part of the spectrum around the pump. The merging of 
the spectral parts around A as , A p , and A s would create 
an ultra broad spectrum as observed in tapered fibers 
with high power femtosecond pulses [15,16]. However, 
for high power femtosecond pulses the SPM is the dom- 
inant mechanism that leads to broadening and merging 
of the Stokes and anti-Stokes lines with the pump. For 
low power picosecond pulses the Raman and parametric 
processes are dominant. Thus, in this particular case, 
the large frequency shift ~ 193THz of the degenerate 
FWM sidebands and the narrow degenerate FWM gain 
bands ~ 2THz prevent merging of the pump with the 
Stokes and anti-Stokes lines to happen within the max- 
imum length L max , i.e., before nonlinear effects become 
negligible. Indeed it is seen from Fig.l, that the spec- 
trum does not change significantly from lm to 3.7m. The 
power transferred to the Stokes and anti-Stokes lines is 
in effect lost, i.e., the SC process is not very efficient. 
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FIG. 2. Dispersion profile dl (thick solid line), d2 (dashed 
line), d3 (dash-dotted line), d4 (dotted line), d5 (circles), and 
d6 (thin solid line). 



By increasing the peak power, it would be possible to 
achieve merging of the spectral parts around A as , X p , 
and A s . However, our aim is to keep the low peak power 
fixed and instead achieve this merging only by engineer- 
ing the dispersion profile. Thus we modify the dispersion 
profile to adjust Eq. (3) to be fulfilled for wavelengths 
A s and A as closer to the pump A p = 647nm. We do 
so by modifying /3 2 , 04, and (3q as listed in Table I. The 
phase-matching condition A/3 = then gives A s = 852nm 



and A q 



521nm for case d2-d3. In case d4-d6 addi- 



tional Stokes and anti-Stokes waves exist. The dispersion 
profiles and phase-mismatch curves corresponding to the 
cases in Table I are shown in Fig. 2 and Fig. 3, respec- 
tively. It is important to note that the curves A/3(A) have 
different slope around A s and A as (see Fig. 3). 
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We first consider only a shift of the Stokes and anti- 
Stokes lines closer to the pump and the effect it has on 
the improvement of the SC generation. Thus, for the 
dispersion profile d2, the slope of the phase mismatch 
curve around A s = 852nm and A as = 521nm is kept the 
same as for case dl. It is seen from Fig. 1 (right) that 
such a shift of the direct degenerate FWM Stokes and 
anti-Stokes lines closer to the pump is not enough for a 
complete merging to take place. This can be explained by 
considering the direct degenerate FWM gain g(X) shown 
in Fig. 4. The broadening of the Stokes and anti-Stokes 
lines is strongly influenced by the bandwidth of g(X), 
which is mainly determined by the slope around A/3 = 0, 
i.e., around A s and A as . The slope and thus the gain 
bandwidth is the same in case dl and d2, which explains 
why the broadening appears to be unchanged. 
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FIG. 3. Phase mismatch A/3(A) for dispersion profile dl 
(thick solid line), d2 (dashed line), d3 (dash-dotted line), d4 
(dotted line), d5 (circles), and d6 (thin solid line). 

One way to improve the SC spectrum is to shift the 
Stokes and anti-Stokes lines even closer to the pump, 
keeping the slope of the phase mismatch curve around 
them constant. However, this will not significantly im- 
prove the width of the SC spectrum as compared to 
case dl, because the narrow direct degenerate FWM gain 
bands will then be in the region, where a SC is already 
generated by Raman and FWM processes. Moreover, 
this will lead to even more unusual dispersion profiles 
than that for case d2 (see Fig. 2), which does not seem 
to be experimentally realistic. 

Instead we fix A s and X as while reducing the slope and 
thus increasing the gain bandwidth. For dispersion pro- 
file d3 the direct degenerate FWM gain bandwidth is 
thus increased to 16.5THz (see Fig. 4). This leads to 
broader Stokes and anti-Stokes lines in the initial stages 
of the SC generation and finally to their merging with the 
spectrum around the pump as seen from Fig. 5 (right). 



Thus a SC, which is flat within lOdB and extending from 
around 500nm to 900nm is formed after a propagation 
distance of L = 2m despite using low-power picosecond 
pulses. Moreover, the dispersion profile d3 is more real- 
istic, i.e., closer to a real PCF dispersion profile, such as 
profile dl. 




(dl) and (d2) 



0- 



400 



1000 



1200 



600 800 
Wavelength [ nm 1 

FIG. 4. Degenerate FWM gain g(X) for dispersion profile 
dl (thick solid line), d2 (dashed line), d3 (dash-dotted line), 
d4 (dotted line), d5 (circles), and d6 (thin solid line). 
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FIG. 5. Dispersion profiles d2 (left) and d3 (right). Top 
row: phase mismatch A/3 for direct degenerate FWM (dashed 
line) and the spectrum at L =10cm (solid line). Bottom row: 
spectrum at L=20cm, 30cm, lm, and 2m (down to up). 

The SC process is thus much more efficient with dis- 
persion profile d3 than with dl, since the power in the 
Stokes and anti-Stokes lines is not lost. However, the 
SC may be further improved by designing the dispersion 
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such that the phase-mismatch A/3(A) has two or even 
three sets of Stokes and anti-Stokes lines, i.e. roots of the 
polynomium (3). The dispersion profiles d4 and d5 rep- 
resents such cases with two and three sets of Stokes and 
anti-Stokes lines, respectively, around which the spec- 
trum can broaden. From the corresponding gain curves 
in Fig. 4 we see that two gain bands actually overlap 
and form one broad gain band. The presence of extra 
Stokes and anti-Stokes lines and the broad gain band 
could make the SC generation more efficient, provided 
they do not deplete the pump so much that the central 
SC spectrum aroud the pump deteriorates. 

From Fig. 6 (left) we see that with the dispersion pro- 
file d4 the initial stage of the SC generation is indeed im- 
proved, the spectrum at L=10cm mainly reflecting the 
gain profile seen in Fig. 4. However, the small dip in 
the gain curve around 950nm has a strong effect on the 
evolution of the spectrum, leaving a clear dip at 930nm 
in the final SC spectrum. Optimizing the position of the 
two Stokes lines can remove this dip and lead to an even 
broader SC spectrum than observed in Fig. 5 for one set 
of Stokes and anti-Stokes lines. 
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FIG. 6. Dispersion profiles d4 (left) and d5 (right). Top 
row: phase mismatch A/3 for direct degenerate FWM (dashed 
line) and the spectrum at L=10cm (solid line). Bottom row: 
spectrum at L=20cm, 30cm, lm, and 2m (down to up). 

Instead we show in Fig. 6 (right) the evolution of the 
spectrum in a PCF with the dispersion profile d5, which 
has three sets of Stokes and anti-Stokes lines. The small 
dip in the gain curve around 800nm (see Fig. 4) is still 
reflected in the spectrum, but it is now less pronounced 
and we obtain a final ultra-broad SC spectrum ranging 
from 450nm to 1250nm within the 20dB level. Of course 
the dispersion profile may be optimized further to remove 
the dip and make the SC spectrum more flat and even 
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FIG. 7. Dispersion profiles d3 (left) and d6 (right). Top 
row: phase mismatch A/3 for direct degenerate FWM (dashed 
line) and the spectrum at L=10cm (solid line). Bottom row: 
spectrum at L=20cm, 30cm, lm, and 2m (down to up). 

So far we have only considered the Stokes and anti- 
Stokes lines generated directly from the pump through 
degenerate FWM. However, so-called cascaded FWM 
processes play also an important role in the evolution of 
the spectrum, as discussed thoroughly in [6] . In particu- 
lar, one can use these processes to obtain a broader SC. 
The dispersion profile d6 is designed to clearly illustrate 
this effect. It still implies two sets of Stokes and anti- 
Stokes lines, but now these are very close to the pump, 
within the regime of wavelengths covered by the central 
SC generated by the pump. Nevertheless we see in Fig. 7 
(right) that additional lines are generated, around which 
the spectrum broadens, resulting in a final SC spectrum 
extending from 450nm to Ifim within lOdB. The line at 
1030nm is generated by direct degenerate FWM from the 
Stokes wave at 720nm, and the line at 1060nm is gener- 
ated by FWM between the Stokes wave at 720nm and the 
pump. Thus these cascaded parametric processes result 
in a spectrum, which is broader than what was obtained 
with the direct degenerate FWM process in case d3. 

Many investigations on designing the dispersion pro- 
file in PCFs have been made [11-13]. In [12] a well- 
defined procedure to design specific predetermined dis- 
persion profiles is established and it is shown that it is 
possible to obtain flattened dispersion profiles giving nor- 
mal, anomalous, and zero dispersion in both the telecom- 
munication window (around 1.55/xm) and the Ti — Za 
laser wavelength range (around 0.8[im). This allows us 
to conclude that the dispersion profiles d3-d5 shown in 
Fig. 2 could indeed be fabricated. However, it is out- 
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side the scope of this work to consider how the specific 
dispersion profiles may be fabricated. 



IV. ROBUSTNESS OF SCG TO FIBER 
IRREGULARITIES 

The SC generation process that we have considered 
here is mainly determined by parametric FWM, which re- 
quires phase-matching. In the experiments with a PCF 
with dispersion profile dl [6,7] Stokes and anti-Stokes 
lines due to direct degenerate FWM were generally not 
observed. This was explained to be due to irregularities 
along the PCF, leading to violation of the required phase 
matching condition (3). Thus, in order to conclude that 
parametric FWM can be used for broad-band SC genera- 
tion in real PCFs, we check the robustness of the process 
towards fluctuations of the dispersion coefficients along 
the fiber. 

It has indeed been shown experimentally that a vari- 
ation of the zero-dispersion wavelength X z of the order 
of O.lnm over the entire length of a dispersion shifted 
fiber could significantly reduce the FWM efficiency [20]. 
This reduction in the FWM efficiency was later explained 
theoretically from expressions for the average parametric 
gain, phase-conjugation efficiency, and gain band-width 
[21]. It has also been shown that in order to control the 
dispersion within ±lps/ (km-nm), the allowable deviation 
of the core radius in W-type dispersion-flattened fibers 
is 0.04/zm, while for other types of dispersion-flattened 
fibers the allowable core radius deviation is 0.1/im [22]. 
As PCFs have an even more complex structure, strong 
fluctuation of the fiber dispersion could be expected too. 
However, to our knowledge a thorough study of the in- 
fluence of fluctuations of the PCF parameters (e.g., core 
size and pitch) on the variation of the dispersion profile 
(i.e. the dispersion coefficients ^2-7) is not available. 

For a newly developed highly nonlinear PCF with 
A z = 1.55/im, it was recently shown that the variation of 
X z is only 6nm and the variation of the dispersion slope 
at X z , D' z = dD{X z )/dX, varies between -0.25 and -0.27 
ps/(km-nm 2 ) over a 150m span [13]. Expanding the prop- 
agation constant to third order around the pump wave- 
length A p =647nm, the dispersion has the form 



h 2tt< ( \ y 
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,/From this expression we find the dispersion coefficients 

XtDl ( 1 1 \ „ AfiX 
X 
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which gives the extrema /3 2 ™ aa; =7.51ps /km, /3™ m = 
6.83ps 2 /km, /?;f a;l; =0.44ps 3 /km, and /?£" n =0.40ps 3 /km 
and thus the relative variations (fo) — (/?3)=9.5%, where 



max 

fc 



P. 



ruin 
k 



{(3f ax +l3f m )/2' 



k = 1,2. 



(7) 



Note that the relative variations of (3i and (3% are equal. 

We model the effect of a fluctuating dispersion profile 
by imposing that 5(3, A, and all the dispersion coefficients 
/?2-7 [see Eq. (1)] vary randomly along the fiber, 

5(3(z) =5(3 + <r (z), Pk,x(z) = (3 k + <Jk,x{z), 
A(z) = A + 0-1(2:), f3 k , y (z) = f3 k +<Tk, y (z), 

where k=2,..,7. The random fluctuation of the coeffi- 
cients, represented by the a-terms, is modelled as Gaus- 
sian distributed white noise with zero mean. To achieve 
the most severe case we use different seeds for all terms. 
We have thus assumed that the fluctuations affect the dis- 
persion in the two bircfringent axis independently. With 
the results from Ref. [13] in mind we assume that the 
strength (or variance) of the fluctuations is the same in 
all coefficients, 



{a 2 (z)) {a 2 (z)) _ _ «„(*)) 



5(3 



(3k 



(8) 



and use the strength p=10%. 

Random fluctuations of the whole dispersion profile 
will randomly vary not only the zero-dispersion wave- 
length X z , but more importantly, the phase-mismatch 
curve A(3(X), given by Eq. (3). This in turn implies that 
the FWM gain spectrum g(A(3), given by Eq. (4), will 
vary randomly along the fiber, even in the undcplctcd 
pump approximation (constant peak power I p ). 

In Fig. 8 we have depicted the fluctuation of the FWM 
Stokes gain band in the undeplctcd pump approximation 
over the first L=lmm of a PCF with dispersion profiles 
dl (A s =1110nm) and d3 (A s =852nm). As expected the 
broader gain band of fiber d3 is reflected in a suppression 
of the oscillations, as compared to the fiber dl used in 
the experiments in [6]. 
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FIG. 8. Random fluctuations of the FWM Stokes gain band 
<?(A/3) (grey region), given by Eq. (4), for constant pump I p 
and p=l% (left) and p=W% (right). The upper row is for 
case dl and the bottom row for case d2. 



The variation of the direct degenerate FWM Stokes 
gain band gives an impression of the influence of the dis- 
persion fluctuations on the effectiveness of the FWM pro- 
cess. The important factor is the average FWM gain over 
the fiber length L, defined as 



9av = j J 5[A/3(z)] dz, 



(9) 



where we have indicated the z-dependence of the phase- 
mismatch as a result of the fluctuations. In Fig. 9 we 
show the average FWM Stokes gain calculated over the 
first L—2cm using the undepleted pump approximation. 
The reduction of the average gain for increasing strength 
of the fluctuations can be clearly observed. 

Theoretically we thus predict that in fiber dl realistic 
fluctuations would significantly suppress the Stokes and 
anti-Stokes lines generated from the pump by direct de- 
generate FWM, as also stated by Coen et al. [6]. The cor- 
responding simulations, presented in Fig. 10 (left), con- 
firms this prediction. Using direct degenerate FWM to 
generate an ultra-broad SC in the particular fiber dl is 
therefore not realistic. 

In contrast, with our proposed fiber d3 with a broad 
gain band, even fluctuations with p=W% should not sig- 
nificantly reduce the FWM effectiveness, which is also 
confirmed by our simulations shown in Fig. 10 (right). 
In our proposed fibers d4-d5 the FWM gain band is even 
broader, indicating that fluctuations will have even less 
impact. Thus our numerical results show that using di- 
rect degenerate FWM to generate an ultra-broad SC in 
really a realistic option. 
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FIG. 9. Average FWM Stokes gain g av over L=2cm, as 
given by Eq. (9) in the undepleted pump approximation. 
Shown is case dl (left) and d3 (right) for p=0 (solid), p=l% 
(dashed), and p=W% (dotted). 
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FIG. 10. Dispersion profiles dl (left) and d3 (right). Spec- 
trum under influence of fluctuations with strength p=0 (top 
row) and p=10% (bottom row) at L=20cm, 30cm, lm, 2m, 
and 3.7m (down to up). 



V. CONCLUSION 

We have numerically investigated SC generation in 
birefringent PCFs using 30-picosecond pulses with 400- 
kilowatt peak power. Our results show that by properly 
designing the dispersion profile and using the generation, 
broadening, and final merging of widely separated pump 
and FWM Stokes and anti-Stokes lines the SC gener- 
ation efficiency can be significantly improved, resulting 
in a broader SC spectrum and a reduced loss of power 
to frequencies outside the SC. Thus, by optimising the 
dispersion profile, we have generated an ultra-broad SC 
ranging from 450nm to 1250nm within the 20dB level. 

We have shown that the key issue is to make sure 
that the Stokes and anti-Stokes lines are generated close 
enough to the pump to be able to broaden and merge 
with the central (pump) part of the SC before nonlinear 
processes are suppressed due to fiber loss and temporal 
walk-off. We have also shown that this in turn requires 
the FWM gain band to be sufficiently broad, which is an 
essential property of our designed dispersion profiles. 

We have finally investigated the robustness of the SC 
generation process in our proposed fibers towards fluc- 
tuations in the parameters along the fiber. Such fluctu- 
ation could be detrimental to the phase-sensitive FWM 
process, which depends on the degree of phase-matching. 
Simulations including random fluctuations of the disper- 
sion profile along the fiber show that the broad FWM 
gain band of our proposed fibers improve the robustness 
and that the process of efficient SC generation survives 
random fluctuations of a realistic strength. 
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